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We study the thermodynamic properties of a partiaUy 
ionized hydrogen plasma in strong magnetic fields, B ~ 
lO^'^ — 10^^ G, typical of neutron stars. The properties of 
the plasma depend significantly on the quantum-mechanical 
sizes and binding energies of the atoms, which are strongly 
modified by thermal motion across the field. We use new fit- 
ting formulas for the atomic binding energies and sizes, based 
on accurate numerical calculations and valid for any state of 
motion of the atom. In particular, we take into account de- 
centered atomic states, neglected in previous studies of ther- 
modynamics of magnetized plasmas. We also employ analytic 
fits for the thermodynamic functions of nonideal fully ionized 
electron- ion Coulomb plasmas. This enables us to construct 
an analytic model of the free energy. An ionization equi- 
librium equation is derived, taking into account the strong 
magnetic field effects and the nonideality effects. This equa- 
tion is solved by an iteration technique. Ionization degrees, 
occupancies, and the equation of state are calculated. 

PACS numbers: 52.25.Kn, 05.70.Ce, 95.30. Qd, 97.60. Jd 



I. INTRODUCTION 

Magnetic fields B - 10^^ - 10^^ G typical of iso- 
lated neutron stars qualitatively modify many physical 
properties of matter (see Refs. [0j2) for reviews). In 
this paper we calculate the thermodynamic properties 
of a strongly magnetized hydrogen plasma at temper- 
ature T ^ 10^ ° — 10^'^ K, which may compose outer 
neutron-star envelopes . As we shall see, the plasma 
under these conditions can be partially ionized, and 
the quantum-mechanical properties of both free elec- 
trons and bound species (primarily hydrogen atoms) are 
strongly modified by the field, which thereby affects the 
thermodynamics. 

The motion of charged particles in a magnetic field is 
quantized into Landau orbitals. The magnetic field is 
called strongly quantizing if the free electrons populate 
mostly the ground Landau level 1^ . This is the situation 
which we are especially interested in. It occurs when the 
electron cyclotron energy hujc = TieB / (mec) (where fi, e, 
me and c are the Planck constant, electron charge, elec- 
tron mass and speed of light, respectively) exceeds both 



the thermal energy ksT and the electron Fermi energy 
ep — that is for temperatures T <ti Tb and densities 
p < Pb, where 



Tb = 3.16 X 10-^ 7 K, pB ^ 0.809 7^/^ g cm" 



(1) 



(see Sec. IIIAl). Here, the parameter 7 — 
h^B/{mlce^) = B/(2.35 x 10^ G) is the electron cy- 
clotron energy in atomic units. 

We will refer to a strong magnetic field when 7^1. 
A number of studies of the equation of state (EOS) of 
matter in strong magnetic fields were based on various 
modifications of the Thomas- Fermi approximation 
This approximation works reasonably well at large p and 
for large ion charge Zi. Abrahams and Shapiro B cs- 

— 1/2 

timate its validity range as p ^ pbZ- ' . We consider 
Zi = 1 and lower densities, for which atoms are present 
in the plasma and contribute to the EOS. 

The atom in a strong magnetic field 7 S> 1 is com- 
pressed in the transverse directions to the size of the 
"magnetic length" : 



(2) 
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where oq = / {mf.e'^) is the Bohr radius. The ground- 
state binding energy grows logarithmically with B and 
exceeds the ground-state energy of the field-free atom by 
order of magnitude at S ~ 10^^ G Q. Ionization equilib- 
rium of atoms in strong magnetic fields was first discussed 
by Gnedin et al. [|o| and Khersonskii ||ll|. Kherson- 
skii |l^ considered also dissociation equilibrium of Il2^ 
species. However, these pioneering works neglected mod- 
ifications of the atomic properties caused by the thermal 
motion of the atoms across the field. 

The motional modifications of quantum-mechanical 
characteristics of the atom arise from the coupling be- 
tween the center-of-mass motion across the field and the 
relative electron-proton motion 113 -1^]. The role of these 
effects was appreciated by Ventura et al. who, how- 
ever, did not treat them quantitatively. An increase of 
the nonionized fraction caused by the motion effects was 
mentioned by Pavlov and Meszaros ||l^, who used per- 
turbation theory applicable to atoms only slightly dis- 
torted from their rest-state cylindrical shape. Quantum- 
mechanical calculations of binding energies and wave 
functions of hydrogen atoms in any states of motion in 
the strong magnetic fields have been carried out only re- 
cently |l|,0. 

Lai and Salpeter |19 2^ evaluated the effects of mo- 



tion on the ionization equilibrium using an approxima- 
tion for the binding energies of moving atoms which does 
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not apply to the so-called decentered states, for which 
the electron-proton separation is large ||l^ , |l5| , |l7|] . Non- 
ideality effects were included in the ionization equilib- 
rium equation only as a pressure-ionization factor for 
p ^ 10^ gcm^'^. As a result, this equation contains a 
factor which diverges (and becomes even negative) at suf- 
ficiently high temperatures. 

Recently, Steinberg et al. calculated the second 
virial coefhcient of the proton-electron plasma in arbi- 
trary magnetic field and constructed an EOS at low 
densities. The bound states were included using the 
Planck-Larkin partition function. This approach yields 
correct EOS at the low density where the virial expan- 
sion holds |2^. However, the Planck-Larkin formalism 
fails at higher densities, where one has to resort to the 
chemical picture of the plasma, as discussed in detail, 
e.g., by Dappen et al. [p3| . In addition, atomic binding 
energies were calculated in pH using approximations ||Tg| ] 
which have very restricted applicability as shown in p4[ . 

In this paper we use new fitting formulas to atomic en- 
ergies and sizes [ p4[ based on a previous numerical study 
[ p^t , valid for any state of atomic motion. The molec- 
ular II2 fraction is evaluated following the approach of 
Lai and Salpeter but with a modified treatment of 
nonideality. Our knowledge of the quantum-mechanical 
properties of molecules in a strong magnetic field is still 
incomplete, but the evaluation of the molecular fraction 
is useful to determine the validity domain of our EOS at 
relatively low temperatures (where the molecules domi- 
nate). 

The next section presents a simple thermodynamic 
model of the hydrogen plasma. The model is tested in 
the nonmagnetic case by comparison with more elaborate 
models, and is shown to provide sufficient accuracy at 



high T where the molecular fraction is small. In Sec. HI 



we consider a fully ionized plasma in a strong magnetic 
field. The partial ionization and dissociation are dis- 
cussed in Sec. pv| , where an analytic model of the plasma 
free energy is constructed and the ionization equilibrium 
equation is derived. Numerical results are presented and 
discussed in Sec. 0. 



II. THERMODYNAMIC MODEL: THE 
ZERO-FIELD CASE 

A. Chemical picture of the plasma 

A theoretical description of partially ionized plasmas 
can be based either on the physical picture or on the 
chemical picture of the plasma |23| . In the chemical pic- 
ture, bound species (atoms, etc.) are treated as elemen- 
tary entities along with free electrons and nuclei. In the 
physical picture, nuclei and electrons (free and bound) 
are the only fundamental constituents of the thermo- 
dynamic ensemble. The relative merits of the two ap- 
proaches have been discussed, e.g., in [p5|,p6|. 



We use the so-called occupation probability formalism 
in frames of the chemical picture. Occupation probabil- 
ities, which ensure convergence of the internal partition 
functions (IPF), were first introduced by Fermi ||27| , who 
has demonstrated their immanent relation to a nonideal 
contribution in the Helmholtz free energy. Various ap- 
proaches to the construction of the occupation probabil- 
ities have been reviewed by Hummer and Mihalas [E8|. 
The approach adopted by Mihalas and co-workers pS^pG] 
(hereafter MDH) is based on the Inglis- Teller criterion of 
Stark broadening conventional for plasma spectroscopy, 
which gives optical spectra consistent with available ex- 
periments (see, e.g., Ref. [Q). However, the equation of 
state derived by MDH is unrealistic at p 10^^ g cm^^ 
(see H]), and the approximations made in its derivation 
are lacking in self-consistency An alternative EOS 
was derived in a self-consistent manner by Saumon and 
Chabrier [ p6pl|j3^ , p^ (hereafter SC) from effective pair 
potentials between plasma particles, but with neglect of 
the Stark broadening. The ionization degree deduced 
by SC strongly differs from that by MDH. The origin of 
the discrepancy is rooted in the fact that strongly per- 
turbed atoms, whose spectral lines disappear due to the 
Stark merging, may still contribute to the EOS as bound 
species [|5). Thus the approaches of MDH and SC are 
reconciled by an approximate treatment of the atoms per- 
turbed by plasma ions as quasicontinuum atomic states, 
which contribute to the EOS as atoms although they do 
not show atomic spectral lines p^ . 

The chemical picture faces a principal difficulty in cases 
where the interaction between nuclei and electrons in a 
bound state is comparable to the interaction between a 
bound object and neighboring plasma particles. This 
situation occurs when pressure ionization is important 
or when high atomic levels are appreciably populated. 
In these special term should be included into 

the free-energy model, in order to distinguish between 
bound and free states. For instance, MDH constructed 
an ad hoc "pressure ionization term" in the free energy 
pof , SC introduced hard cores with fixed diameters in 
the effective potentials for bound species js^, and expo- 
nential "unbinding" occupation probabilities were used 
in 1^^. The latter approach has been justified by consid- 
ering an excluded volume of the bound objects at rela- 
tively low density, assuming an uncorrelated distribution 
of the plasma particles. At high density, the strong cor- 
relations of the positions of the particles must be taken 
into account. Their approximate treatment in the hard- 
sphere model (e.g., by SC) appears to be practical for 
this purpose. 

In the case of the strong magnetic field, the model of 
the plasma cannot be as detailed as, e.g., the SC non- 
magnetic model, because the effective potentials (partly 
derived from high-pressure experiments in the zero-field 
case) are not available. Therefore we use a simple hard- 
sphere picture described below. In order to check the 
validity of this model, we apply the same assumptions to 
the well-studied zero-field case and compare the results 



2 



with those of more elaborate models. 



B. Free-energy model 

Consider a plasma consisting of electrons, protons, and 
H atoms in a volume V. Let us write the Helmholtz free 
energy as = Fjd + Fcx, where 



id 



F: 



Fmcu 
id 



Fr. 



ad 



(3) 



is the sum of the ideal-gas free energies of the electrons, 
protons, neutral species, and photons (thermal radia- 
tion), respectively, and F^x is the excess (nonideal) part. 



1. Ideal part of the free energy 

We consider nondegenerate protons and neglect their 
spin statistics both in bound and free states (this 
amounts to an additive constant in the entropy that af- 
fects neither ionization equilibrium nor the EOS, pro- 
vided the total number of free and bound protons iVo is 
fixed). Then 



(4) 



where f3 = [kBT)"^. Here and hereafter, TV^, n^, and 
Aa = (27r/9?i^/mQ)^/^ denote, respectively, the total num- 
ber, number density, and thermal wavelength of particles 
of the type a with mass rria . 

For the ideal gas of electrons, we use the identity pq] 



F: 



id 



(5) 



where and Pg are the chemical potential and pres- 
sure of the ideal Fermi gas, respectively, which can be 
obtained as functions of Ue and T from the equations 



P. 



8 ksT 



4/2 (/^Me), 



A3 

^/l/2(/3/.e). 



Here, 



exp(a; — z) -|- 1 



da; 



(6) 
(7) 

(8) 



is the Fermi integral. With the use of Fade approximants 
to the functions and their inverse functions [ p7| , 

-^id^ is expressed as an analytic function of A^g, V, and 
T. 

In the zero-temperature limit, one may replace I^{Pfie) 
by (/3ei?)'^+^/(i/ -I- 1), which gives, in particular, the well- 
known expression 



2mf 



(37r2n,)2/3. 



(9) 



The Fermi temperature is defined as Tp = ei^/fes « 3 x 
105p2/3 where p = 1.6735 ne/(10^'' cm^^-j ^j^g j^^ss 
density of the electron-proton plasma in gcm"^. in the 
nondegenerate limit T ^ Tp, the ideal Boltzmann gas 
relations are recovered, fie — fcBnn(neAg/2) and Pe = 
neksT. 

For the atoms, one has 



= ksTj^N. [ln(n,A3/g«) - 1 



(10) 



where n enumerates quantum states with statistical 
weights and binding energies Xk - 

It should be noted that, although nonideality effects 
are not included in Fid explicitly, they do affect the equi- 
librium value of Fid through the particle numbers. In 
particular, the distribution of A^^ in Eq. ( p^ is not as- 
sumed to obey the ideal-gas Boltzmann law. 

Finally, the radiation term (which can be important 
only at low p or very high T) reads 



-Frad = -(4a/3c) VT\ 



(11) 



where a = t: kg/{60h c ) is the Stefan-Boltzmann con- 
stant. 



2. Excess free energy 



The excess free energy is conventionally written as 



(12) 



where F^ is the excess free energy of the ionized part of 
the plasma and F™" accounts for interactions of neutral 
species with electrons, protons, and other neutral species. 
The Coulomb term 



7C 



F, 



(13) 



includes contributions from the exchange and correla- 
tion interactions of electrons Fee, Coulomb interactions 
in the one-component plasma (OCP) of ions Fa, and ion- 
electron (screening) interaction F^g. These contributions 
have been calculated by various procedures, e.g., by solv- 
ing a set of hypernetted-chain equations or A^Ionte Carlo 
simulations pSj-fl^ . We make use of the fitting formulas 
to the results of such calculations, obtained in for Fee 
and in for Fa and Fie- These formulas express the 
electron-ion plasma free energy as an analytic function 
of the electron density parameter 



tte/ao « 1.39p 



-1/3 



and Coulomb coupling parameter 

r = /3eVae« 0.227 pi/ VFe, 



(14) 



(15) 



where Oe — (47rne/3) ^/"^ jg the mean interelectron dis- 
tance and Tg = K. 
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The nonideal part of the atomic free energy, F^- 
be written as |33y34| ] 



neu 

cx ' 



pnCU 771 



pert ; 



(16) 



where Fus is the reference free energy, treated in the 
hard-sphere approximation, and -Fpert is the perturbation 
part that accounts for the attractive (van der Waals) in- 
teractions. To calculate these contributions, an elaborate 
model has been developed by SC [ ^J34| . Its simphfied 
analytic version for weak electron degeneracy has been 
constructed in |Q. In the so-called van der Waals one- 
fluid model 1^1, a free energy of the hard-sphere mixture 
is represented by the Carnahan-Starling formula |44| 



l3F^s/Ntot = i^V-^V^)/{^-vf, 



(17) 



where A'tot = J2a is the total number of particles, 



GNtotV 



(18) 



is the effective packing fraction, and daa' are the hard- 
sphere interaction diameters. In our case, the subscript 
a enumerates atomic quantum states described by quan- 
tum numbers k and takes on a single value p for the free 
protons. 

In the following, we compare two versions of the model: 
(i) the full version, in which -Fport and daa' are given 
by ap pro ximations of with one exception (adopted 
from [p4| ) that daa' cannot be smaller than a certain 
limit d''^^,, and (ii) the simple version, in which long- 
range atomic interactions are disregarded. In the latter 



pert 



and da 



7(0) 



Furthermore, we adopt 



the simplest choice 



- I 



- I 



(19) 



where 1^ is the root-mean-square proton-electron dis- 
tance in the quantum state n p5[ . For the interactions 
among charged particles, we define daa' ~ 0, because this 
type of interaction is already included in the Coulomb 
part of the free energy. Note that in the second (simple) 
version of the model, .Fhs turns into the unbinding term 
Fub of Ref. Q in the low-density limit {j] < 1). Thus 
the unbinding term is now incorporated in -Fhsi which 
allows us to approximately take into account the corre- 
lation effects. 



C. Equilibrium conditions 

Thermodynamic equilibrium is given by the minimiza- 
tion of F{V,T, {Na}) with respect to the particle num- 
bers Na under stoichiometric constraints. The condition 
of the extremum of F can be written in the form of the 
Saha equation corrected for nonideality and electron de- 
generacy: 



riH = npneXl{mp/miif^^{Z^/2)e^, 

where 

A = l3dF^^^/dN, - Hn,Xl/2) 
allows for electron degeneracy, and 



(20) 



(21) 



(22) 



is the modified IPF which includes the occupation prob- 
abilities Wk. , defined according to [p2l : 



/cbTIu w„ 



dF,^ dF,, 



dN„ 



(23) 



To solve Eq. ( ]20| ) , one must add the electroneutrality con- 
dition Tie — Up and the mass conservation condition njj -t- 
np — uq, where no = p/nin = (p/11.293 gcm^"^) Uq^. 

The Boltzmann distribution of the atoms, corrected for 
nonideality, reads 



0x. 



(24) 



The minimum of the free energy is sought by solving 
Eqs. (20)~(|2^) iteratively First, one defines star ting 
Wk's and calculates the number densities from Eqs. ( |2(]| ) 
and (|2^). Then the w^'s are refined using these number 
densities in Eq. (H) @. 

The molecules II2 can be easily included in this proce- 
dure. The dissociation-recombination equation reads 



?1H: 



, =<(Ah\/2)3Z^2A 



2 

w 7 



(25) 



where is the molecular IPF, modified by multiplying 
each Kth term by an occupation probability [^ , 
given by 



kBTlnw^^ = 2 



fdF,., 



\dNp 



dF^ 



(26) 



For simplicity, we do not include molecules in the present 
versions of the model, because the fraction of II2 is small 
in the range of p and T which we are interested in. 

After the equilibrium distribution of plasma particles 
is found, the pressure P, internal energy J7, and entropy 
S are calculated from the relations 



P 



-(aF/9y)T,{jv„}, u ^[dim/dp] (27) 



S — (U — F)/T. The higher-order thermodynamic quan- 
tities are obtained by differentiation of P, U, S without 
keeping Na fixed 



D. Results of comparison 

The ionization curves given by different versions of 
the model are compared in Fig. |l] for T = 10^'^ K. Al- 
though the neglect of the perturbation terms introduced 



4 



in the simple version is most perceptible at such rela- 
tively low temperatures, the "full" and "simple" versions 
of the model yield practically identical atomic fractions 
/h = nu/nn. 




-3 -2 -1 

logloP (g cm-=) 

FIG. 1. Comparison of the fraction of neutral atoms 
/h = nu/no given by SC tables (short-dashed line) and by 
two present versions of the thermodynamic model of partially 
ionized atomic hydrogen (see text). The long-dashed line cor- 
responds to the fraction of atoms that satisfy the Inglis- Teller 
criterion. The dotted line is given by the usual Saha equation 
for ground-state atoms. 

The results of SC qualitatively agree with the 
present model. Quantitatively, they differ in the 
pressure- ionization region at p > 0.1 gcm~^, where the 
theoretical uncertainty is largest (see Sec. [I A). The 
difference in the low-density regime p < 10"'' gcm^'^ is 
due to highly excited states. If both the effective diame- 
ter d and statistical weight gn are proportional to n'^ {n 
being the principal quantum number), then the neutral 
fraction should asymptotically decrease at low density as 
/h oc p^/^. Our present model exhibits this asymptotic 
behavior; the dependence /h oc p^^^ seen at low p in the 
SC data might result from a choice d oc n. 

The long-dashed curve represents the fraction of atoms 
satisfying the Inglis- Teller criterion: /it = "Erii^WK/nQ. 
Here, Wk is the probability that a given atom is not 
strongly perturbed by plasma microfields; it is estimated 
from Eq. (31) of Ref. Using /it, we have calculated 
monochromatic opacities of the plasma and compared 
them with the OPAL monochromatic opacities (at 
p < 10 r| gcm~^ where the OPAL data exist). Along 
the isotherm shown in Fig. ^ our results agree with 
OPAL within 12% in the photon energy range from 13.6 
eV to 500 eV where the opacity is dominated by bound- 
free atomic absorption. For comparison, hydrogen opaci- 
ties calculated in ^ differ from OPAL by up to 37% (in 
the same range of energy and density at the same T) . 




logioP (g cm-3) 

FIG. 2. Comparison of two present versions of the EOS of 
partially ionized atomic hydrogen (see text) with Refs. ^] 
(SC) and m (P96). The EOS of ideal fully ionized gas is 
also shown. 

Figure ^ demonstrates that the EOS obtained with the 
full and simple versions of our model practically coincide. 
In the region of weak degeneracy, they also coincide with 
the model presented in Moreover, there is a good 

agreement with the SC model Small differences oc- 
cur only in the regions where the SC model predicts an 
appreciable amount of molecules, as explained in ||3^ . 

As is well known, the second-order quantities are more 
sensitive to the details of the thermodynamic model than 
the first-order ones. The adiabatic temperature gradient 



Vad = (a In T/a In P)s 



(28) 



is shown in Fig. ^. There are only tiny differences be- 
tween the full and simple versions. For comparison, we 
also show Vad from other models. In its validity region 
(i.e., at low density), the model approximately agrees 
with the present one. The differences with predictions of 
SC are somewhat larger. In all models, the isotherms 
"wiggle" in the region of consecutive pressure destruc- 
tion of excited atomic states. Such wiggles are absent in 
the OPAL data ji^ , based on the physical picture of the 
plasma and also shown in Fig. |[ Compared to SC, the 
present data tend to be closer to the OPAL data. We 
conclude that the simplifications introduced above are 
acceptable to describe the thermodynamics of atomic hy- 
drogen. In Sec. we generalize the model to the case 
of the strong magnetic field. 
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FIG. 3. Two adiabatic gradient isotherms given by differ- 
ent EOS models in the domain of partially ionized atomic 
hydrogen. 



III. FULLY IONIZED PLASMA IN A STRONG 
MAGNETIC FIELD 

In this section, we describe effects of quantizing mag- 
netic fields on the fully ionized proton-electron plasma. 
We assume that the field B is uniform and directed along 
the z-axis. 



A. Ideal gas 

1. Electrons 

The electron energy in a magnetic field reads pO[ | 

ew(p.) = Nhujc + p'/(2me), (29) 

where Pz is the longitudinal momentum and N ~ 
0, 1, 2 . . . is the Landau quantum number. All levels ex- 
cept the lowest one are double degenerate with respect 
to the spin projection. Strictly speaking, the anomalous 
magnetic moment of an electron leads to a splitting of 
the levels iV > 1 by 0.00116 hiOc, which takes off the dou- 
ble degeneracy. However, this splitting cannot affect the 
thermodynamics at p < ps, where fcsT should be at 
least comparable to htOc for an appreciable population of 
the excited Landau levels. 

The thermodynamic functions of the electron gas in 
the magnetic field are easily derived from the first prin- 
ciples [p6[ . Taking into account the fact that the number 



of quantum states of an electron with fixed discrete quan- 
tum numbers in volume V per longitudinal momentum 
interval Ap^ equals VApz/iiir'^a'^h) |50|, the thermody- 
namic potential Q = —PV can be written as 



VkuT 



oo 

™ N=0 



9n 



In ( 1 e 



/3[Ale-CN(Pz)] 



dPz, 



where is the statistical weight (go — 1 a-nd gi^ — 2 for 
iV > 1). Hence the electron pressure and number density 
are given by the equations 



1 



27r3/2„2 \ 



(30) 
(31) 



where /i^r = p^ — Nhujc- The Helmholtz free energy is 
given by Eq. (|^), where pe can be found by inversion of 
Eq. (0) (e.g., using an algorithm described in pl|| ). 

In the nonquantizing magnetic field Tg ^ T, where 
many Landau levels are populated, the sum over N in 
Eqs. (^o|), (|3l| ) may be approximated by an integral, and 
integration by parts reproduces Eqs. (^), (0). 

In the domain of strong magnetic quantization, T <C 
Tb and p < pB, one may retain only the term = 0. 
In that case, replacing I-i /2(Ppe) by 2\J (iep in Eq. ( |3l| ) 



(by analogy with Sec. II B) yields the Fermi energy 



ep = 



(32) 



By definition, p — pB &t ep = fioJc- Hence p/ pb = 
"in {2'y rlY'^ / ^ from which Eq. (||) follows. A comparison 
of Eqs. (^) and ( p2[ ) reveals that the Fermi energy changes 
by a factor (4/3p/3(p/pB)4/3. Thus the de generacy is 
strongly reduced at p <C ps. 

In the nondegenerate regime T ^ Tp, one has 
luiPpe) ~ exp(/3^e)r(i/ -I- 1); therefore Eqs. @), (js^) 
reduce to = riekBT and 

(3pe = ln(7ieAe/2) - Inw + In(tanhu), (33) 

where u = (3hujc/2 = Tb/(2T). This yields an ex- 
plicit analytic form for Pj^^ In the nonquantizing field. 



Tb ^ T, the last two terms i n Eg . ( p3[ ) cancel out, and 
the classical expression (Sec. [IB ) is recovered. In the 
strongly quantizing regime, p < ps and Tp <^T <^Tb, 
the last term of Eq. (p3) vanishes, which yields 



P^^^ = NekBT [ln(27ra2„Ae?ie) - l] 



(34) 



2. Protons 

The transverse motion of the protons is quantized in 
Landau orbitals with the elementary excitation equal to 
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the proton cyclotron energy hu^ 



{me/'rnp)fujJc- The 



^2 r 



energy spectrum is given by Eq. ( pQ]) when replacing me 
by rup and lOc by LOcp. Unlike the case of electrons, the 
double spin degeneracy of the Landau levels is taken off 
by the abnormal magnetic moment of the proton. 

In our analysis, the protons are always nondegenerate, 
so that by analogy with Eq. (^) we have 

l3F,f/Np = ln(27ra^Apnp) + In [l - e"^'"^-] - 1. (35) 

Here, for sake of brevity, we drop the zero-point energy 
^TiLUcp and the spin energy zLjgphuicp, where gp — 5.585 
is the proton spin gyromagnetic factor . We suppress 
these terms also for atoms and molecules. Taking them 
into account yields an additive contribution to the total 
free energy of the system, equal to 

AF = No {^nuj,p - fcsT ln[2 cosh{/3gphuJcp/4:)] } (36) 

(iVo = Np in the case of full ionization). Since Nq is 
constant, AF affects neither ionization equilibrium nor 
pressure, but it does affect the internal energy and spe- 
cific heat, therefore we take it into account in Sec. 



B. Nonideal Coulomb plasma 

According to the Bohr-van Leeuwen theorem, mag- 
netic field does not affect the thermodynamics of classi- 
cal charged particle systems (see, e.g., Ref. [|2]). Thus 
the classical ionic OCP excess free energy Fii{r) does 
not depend on B at any F. The classical regime for 
the electron-proton plasma corresponds to ^ 1 and 
r 1, where the excess Coulomb free energy is given by 
the Debye-Hfickel formula Fg = -NeksT ^SE^/S. In- 
deed, it is easy to check that this law holds independent 

ofi?i. 

A magnetic field, however, affects quantum-mechanical 
contributions to F^. These effects have been studied 
only in low-temperature or low-density regimes. 

The ground-state exchange energy of the electron 
gas in the strongly quantizing field behaves as 

-2.25(7r3)-i[ln(7r2) - 0.457 H ] e^/ao (per electron), 

compared to — 0.757r^^(97r/4)^/'^r7^e^/ao in the nonmag- 
netic case Thus the exchange energy at T Tp is 
suppressed by a factor « 0.2036 7r^/[ln(7r^) — 0.457]. 
Note that the condition of strong magnetic quantization 
requires 7r^ to be large. 

A general low-density expansion for the free energy 
of a Coulomb plasma in arbitrary magnetic field up to 
order p^^^ has been derived by Cornu [Q. The coeffi- 
cients of this expansion are not available in explicit an- 
alytic form but require numerical evaluation, which has 
not been done yet. In the particular case of the OCP, a 
Wigner-Kirkwood-type expansion in powers of Ti is avail- 
able . The lowest-order term of the latter expansion 
(quantum diffraction term of order h^) has been obtained 
by Alastuey and Jancovici p5|: 



Fdiff = N.kBT 



8r,, 



tanh? 



(37) 



with u defined as in Eq. ( |33| ) . The expression in the 
square brackets in Eq. ( pT)) goes to 1 at u ^ 0, recover- 
ing the well-known zero-field result, and to ^ at u ^ 1, 
refiecting the fact that two of three degrees of freedom for 
the electron motion are frozen out in the strongly quan- 
tizing field. Equation (|37| ) is valid in the low-density 
regime, where rs ^ max(r,r^^). In this regime the 
correction (37) is smaller than the classical OCP correc- 
tions to the Debye-Hiickel formula. In the electron-ion 
plasmas, Fdiff is canceled because of the local neutrality 
relation [s^ . 

In this paper, we are mainly interested in the case 
where (Fr^)^^ « 3.16 Tg ^ 1- In this case, a high- 
temperature expansion |56[ , which can be written as an 
expansion in powers of two small parameters si = Vt/ts 
and S2 = -s/FrJ, is relevant at low densities. The lowest- 
order correction is the Hartree-Fock term oc h^e^. Stein- 
berg et al. |^l| have recently obtained an analytic result 
for this term in a magnetic field: 



PF, 



HF 



SF^ 

"8r, 



(38) 



where the function 



cosh(2u) 
cosh^ u 



tanh u 



arctanh (^y/l — u ^ tanhu) 



Vl — u tanh u 



(39) 



goes to 1 at small u, reproducing the zero-field result [ p6[ , 
and to ln(4w)/M at very large u. 

Steinberg et al. |^ have also calculated the corrections 
cx he'^ (the MontroU-Ward and exchange terms). For the 
electron gas, they can be written in the form 



/3F4 3V^ F5/2 



N, 



16 



(40) 



where fi'^iu) and f^Hu) go to 1 at u ^ 0, reproducing 
the zero- field result and decrease at large u. 

In order to incorporate these results into the analytic 
free-energy model, we employ a simple scaling procedure. 
In the fit for fee{0,T) = jSFee/Ne derived in (39), where 
9 = T/Tp is the degeneracy parameter, we replace the 
zero-field value 6*0 = 2 {9n/'i)-^/^rs/T by 



1 



(41) 



Here, 9„i = 8'y^r^J{9Tr'^r) = 0.166 6*0 7^^4 is the value 
of the degeneracy parameter in the strongly quantizing 
field. 

The scaling (^Tj) reproduces limiting cases: (i) at 
Ts ^ 1, the classical OCP expression is recovered, inde- 
pendently of other parameters; (ii) in the nonquantizing 
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regime jr^ <C 1, we get the nonmagnetic value 9* = 9q; 
(iii) in the strongly quantizing degenerate regime ^r'^ S> 1 
and ^ 1, the correct value of the degeneracy param- 
eter 9* = 9m is recovered; and (iv) in the strongly quan- 
tizing nondegenerate regime jr'^ ^ 1 and 9m ^ 1, the 
fit reproduces Eq. ( ^8|) in its range of validity. 
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FIG. 4. Contribution to the electron-gas nonideal free en- 
ergy /ee = /3Fee/iVe beyond = -r^^VVS at three val- 
ues of the magnetic field parameter 7 and three values of the 
density parameter Vs (indicated). The scaled fit (solid lines; 
see text) is compared with the high-temperature expansion 
up to orders and (dotted and dashed lines, respectively) 
and with the classical OCP (dot-dashed lines). 

Figure^ demonstrates the validity of the adopted mod- 
ification of Fee at r < 1, for three values of r^, for which 
the quantum contributions to could appreciably af- 
fect our results. We plot departures of Fee from the OCP 
Debye-Hiickel function F§^^ , normalized to NeksT and 
divided for convenience by . The dot-dashed line shows 
the classical OCP free energy the dotted line dis- 
plays the correction (^8[), and the dashed line results 
from inclusion of the corrections (|40|). 

The left panel presents the nonmagnetic case. The 
solid line shows the fit to Fee IH] ■ The region of approxi- 
mate coincidence of the fit with the high-T expansion can 
be adopted as the region of validity of the latter. At large 
Ts, it is restricted by the condition S2 ^ 1 (F ^ rj^), 
while at small r^, the condition si ^ 1 (F <C r^) is more 
restrictive. 

The middle and right panels show the modifications 
of Fee at two values of the magnetic field strength. One 
can see that the scaled fit satisfactorily reproduces the 
expansion in the validity range of the latter. Surprisingly, 
although the scaling is based on the lowest-order term 
(^), the terms ( po| ) are also well reproduced. 

For the electron-ion plasma, the screening contribution 
Fie should be taken into account. At B — 0, it has been 
calculated in a number of papers (e.g., Refs. p9|-pl|) and 



fitted by analytic expressions . In a strongly quantiz- 
ing magnetic field, this contribution has been analytically 
evaluated only for a dense plasma at zero temperature us- 
ing the linear response theory . Comparison with the 
analogous zero-field result [5^ shows that the strongly 
quantizing magnetic field increases the screening energy 
at high density by factor 0.8846 7^rg. To our knowledge, 
there were no relevant calculations at arbitrary degener- 
acy. In the regime of low degeneracy and weak Coulomb 
coupling, integral representations of the low-density ex- 
pansion coefficients have been obtained mJ52] . The con- 
tribution of order is given by Eqs. (20)7(21) of Ref. 
pl[ . It is reproduced if to replace by r^/ (/l^)^ in the 
nonmagnetic expression [S^ . Here, can be approxi- 
mated (within 0.5%) as 



0.9 



arctanh [(1 — t) 
2(l-i)"-6 



0.61 



(42) 



where t = tanh(0.4M)/(0.4u). We apply the scaling ; 
fs/ifTY to the formula for Fje(rs,F) given in Q. 



IV. PARTIALLY IONIZED PLASMA 

A. Hydrogen bound species in the strong magnetic 

field 

1. Atoms 

Only a brief summary of the properties of the hydrogen 
atom in a strong magnetic field is given below. See, e.g., 
Ref. [0 for details and references. 

The motion of an atom in a magnetic field B can be 
conveniently described using the pseudomomentum K, 
the quantum-mechanical constant of motion related to 
the average center-of-mass velocity v ~ Vk-E, where E 
is the total energy of the atom. If there were no Coulomb 
attraction, the energy would he E = E^^ + K^/{2mH), 
where 



E]^.. = Nhcuc + (iV + s)hujcp 



(43) 



is the energy of the transverse excitation, N is the elec- 
tron Landau number, s is the z projection of the rela- 
tive proton-electron angular momentum, and / (2mB_) 
is the kinetic energy of motion along the field. The 
Coulomb interaction mixes the Landau orbitals. Nev- 
ertheless, it is convenient to keep the quantum numbers 
N and s for enumerating the quantum states at 7 S> 1. 
Then the energy of the atom can be decomposed as fol- 
lows: 

i?Ar,,(K) - eI^^{K^) + Ei,, + Kl/(2vin), (44) 

where e\j^^{K^^) < is the "longitudinal" energy, and 
the quantum number v enumerates the longitudinal exci- 
tations. At 7 1, the states with ^ and large s are 
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subject to autoionization. Therefore we put N = and 
suppress this quantum number hereafter. The binding 
energy is 



sfiuj 



cp- 



(45) 



In accordance with Sec. Ill A 2[ the zero-point and spin 
terms are subtracted from Eq. (43) and absorbed into 
Eq. (H). Note that Eq. (||) is vahd for the partially 
ionized plasma provided that the IPFs for atoms with 
opposite proton spin projections are identical. It is true 
under the assumption that the autoionization processes 
with proton spin flip may be neglected on the plasma 
relaxation time scale. We adopt this assumption, be- 
cause the plasma under consideration is rather dense and 
nonrelativistic. Otherwise, states with binding energy 
Xsu{K±) < gphuJcp/2 should be excluded from the IPF 
for atoms with the negative proton spin projection. 

At K = Q, the atom is axially symmetric, and its sizes 
transverse to the magnetic field can be approximated 

as ^ ly K. a. 
size is much larger: 



ys+T, while the longitudinal 
~^ flo/ In 7 for the tightly bound 



states {v = 0) and 1^ ~ ao'^ for the hydrogenlike states 
{v > 1). Longitudinal energies of the former states grow 
as oc (In 7)^, whereas the energies of the latter states 
are relatively small, \E^'^\ ~ (e^/ao)(2n^)~'^, where n is 
the integer part of {v + l)/2. 

An atom moving across the field acquires a constant 
dipole moment in the direction opposite to its guid- 
ing center Yc = c(ei?^)~^B x K. When K±^ is small 
enough, the dipole moment is also small, and E'" is in- 
creased by K\/{2m-^^j). Here, m^^ is the so-called effec- 
tive transverse mass, which exceeds mjj and grows with 
field strength. In this case the average transverse ve- 
locity is v± = K±/mj-j,. When K± exceeds some criti- 
cal value Kc ~ lO^h/ao, the atom becomes decentered: 
w_L reaches a maximum and starts to decrease, while the 
electron-proton separation approaches Tc- Thus, for the 
decentered states, the transverse pseudomomentum K±^ 
characterizes electron-proton separation rather than ve- 
locity. 

In the limiting case of K± ^ ^{v -\- ^)'^h/ao, the lon- 



gitudinal energies approach the asymptote -E" 



Note that only the states with s = may remain bound 
if they have such large values of K±. Indeed, since Ell IS 
small for large K± , the binding energy (^) becomes neg- 
ative for s > 1. However, at s = and arbitrarily large 
Kj_ , there still remains an infinite number of truly bound 
states (enumerated by v), as has been strictly proved in 

Since Tc = aQK±/jTi, the decentered states have huge 
sizes at 7 < 1; hence they are expected to be destroyed 
by collisions with surrounding particles in the laboratory 
and in white-dwarf atmospheres js^. In neutron-star 
atmospheres at 7 > lO'^, however, the decentered states 
may be significantly populated, as we shall see below. 

Accurate numerical dependences of the atomic binding 
energies Xsu{K±) and sizes lsu{K±) for 300 < 7 < 10"* 



and any K± have been obtained in ||l^,gj|. In Ref. p4| , 
analytic fits have been constructed for these quantities 
as functions of 7 and K± , as well as for the critical pseu- 



domomentum Kr and transverse mass 



as functions 



of 7, for various s and i^. These fits are more suitable for 
studying the thermodynamics of hot plasmas than pre- 
vious approximations [ p^ , which were accurate only for 
the centered ground state (s = 1/ = and K± < Kc). 



2. Other bound species 

At sufhciently low T or high B, there may exist a 
considerable amount of molecules and ions in the hy- 
drogen plasma. The molecular ion H2~'' has been thor- 
oughly investigated at i? < 10^° G, including the depen- 
dence of binding energies of various electron- vibrational- 
rotational levels on the angle between the ion axis and 
the magnetic field direction |5^]. For stronger fields, only 
the parallel configuration has been considered ]T2| , |60| -|62[ . 
The ions H2^ have negligible abundance in the strong 
field, owing to the small binding energy, compared to the 
atoms and H2 molecules |l^j22]. The same is probably 
true for the H~ ions l20[| . 

H2 molecules have been studied in detail at various 
field strengths ||60| , |6l| , |63| . An interesting result is that 
the ground state is unbound at 0.18 < 7 < 12.3 [Q. Fit- 
ting formulas for the dissociation energies in the parallel 
configuration for 7 > 10^ have been given in | p0| , |6l[ . At 
such fields, the dissociation energy grows cx (In 7)^, ap- 
proximately at the same rate as the atomic ground-state 
energy. The equilibrium internuclear distance decreases 
as l/ln7, being as small as jao a,t B = 10^^ G, again 
roughly proportional to the longitudinal size of the atom. 

Moreover, strong magnetic fields stabilize polymer 
chains H^r aligned with B, as first suggested by Rud- 
erman J64| and later confirmed by Hartree-Fock calcula- 
tions | |60| |. The specific quantum-mechanical properties 
of these species (e.g. their excitation spectra) are poorly 
known. 

Motional effects on the molecules and chains in the 
strong magnetic fields have not been studied. Therefore, 
one cannot construct a reliable EOS in the domain of 
p, T, B where these species are expected to dominate. For 
instance, Lai and Salpeter estimated the effective 
transverse mass of H^r as iV times the atomic effective 
mass, NniQQ, and used it in the dissociation equilibrium 
equation. However, since the heavier molecule has lower 
velocity at a given K, it is exposed to a weaker electric 
field in the comoving frame. Therefore, one could expect 
its energy levels to be less perturbed and its effective 
mass to be closer to the zero-field value, Nmn. 

Because of these uncertainties, we do not include H^r 
in our study but restrict ourselves to the atomic phase. 
Nevertheless, we include ground-state H2 molecules in 
order to determine the validity domain of our results. 
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B. Free-energy model 

Our free-energy model is a straightforward general- 
izati on t o the magnetic case of the model presented in 
Sec. ETbI: 



rad 



cx ' cx 



(46) 



The ideal electron a nd prot on f ree ene rgies F^^^ and 



are derived in Se es. HII A 1 and III A 2 , respectively. Frad 
is given by Eq. The Coulomb part F^ has been 

discussed in Sec. IIIB. Now let us consider the ideal and 



nonideal contributions and F^^°^ brought about by 
the bound species. 

Since the quantum-mechanical characteristics of an 
atom in a strong magnetic field depend in a nontriv- 
ial way on the transverse pseudomomentum K±, the 
distribution of atoms over K± cannot be written in a 
closed form, and only the distribution over remains 
Maxwellian. Let PsuiK±)d^Kj_ be the probability to 
find an atom with given (s, v) in an element d^K± 
near the point Kj^ of the transverse pseudomomentum 
plane. For the Maxwell distribution, we would have 
Psiy{K±) — (27r?i)~^AH exp[— iir^/(2mH)]. In the general 
case, the number of atoms in an element d^K of the 
pseudomomentum space is 



diV(K) - N,, exp 



PsM<±)<i''K, (47) 



where Ngi, — J diVsy(K) is the total number of atoms 
with the specified discrete quantum numbers. The dis- 
tribution NsuPsuiK±) is not given in advance but should 
be calculated self-consistently by minimization of the to- 
tal free energy, including nonideal terms. 

It is convenient to introduce deviations from the 
Maxwell-Boltzmann distribution through the occupation 
probabilities Wsu{Kj_): 



\2Tih J 

where 



(49) 



/ w;..(i^i)e'5x-^(^-)ifxdi^i, (50) 
2nh^ Jo 



(51) 



The number of atoms per unit phase-space cell equals 
[AN{K)/ d^K]i2'Khf/V. Calculation of {U - TS) for 
this distribution gives 



F^ = ksTY^Nsu j {In [nsM'^^'hfps,{Ki_)\ 



I- l3xsu{Ki_)]ps,{Ki_)d?Ki_ 

3 Wsy{K^) 



(52) 



keTY^Nsu j In 



exp(l)Zs^ 



Ps,{K^)A^K^. 



The contribution of molecules should be added to this 
expression. We estimate it taking into account only the 
molecules in their ground state. This is an acceptable ap- 
proximation at i? > 10^^ G, because the energies of differ- 
ent types of molecular excitations are not much smaller 
than the electronic excitations of the atoms (con- 
trary to the zero-field case), so that excited levels can- 
not give a large contribution to the molecular IPF at 
those relatively low temperatures where the molecular 
fraction is large . We also neglect the (unknown) mo- 
tional modification of the molecular spectrum, as noted 



in Sec. IV A 2 



Thus, we include in the ideal free energy of the bound 
species F^^^ the term 



(53) 



where = 2xoo (0) -t- Q2 is the molecular binding en- 
ergy, and Q2 is the dissociation energy fitted as function 
of 7 in 

The nonideal part F^™ is calculated in the hard-sphere 
approximation using Eqs. (p7|)-(p^), where the composite 
atomic number is k = {svK^_) and the obvious general- 
ization of includes j Psu{K±^)d?'Ki^. The effective 
atomic size Z„ — ls^{K^) is given by fitting formulas [ p4| . 
The effective size of the H2 molecule in the ground state 
is estimated as l^^ = [2a^ + 1^^^]^^'^, where the lon- 
gitudinal size is lz,ii2 ~ ha + fo- Here, Izo is the lon- 
gitudinal size of the ground-state H atom fitted in , 
and tq « 12.7 (In 7)"^'^ is the equilibrium internuclear 
separation given in [po|. 



C. Equilibrium conditions 

The thermodynamic equilibrium for the free-energy 
model ([46|) is given by a generalization of the equations 
in Sec. |ll Q , taking into account the fact that the atomic 
IPF Zu, now includes integration over K±. In partic- 
ular, dFcx/dN^ in Eq. ( P3[ ) is replaced by a functional 
derivative. 

In the conditions studied here, neutral atoms can exist 
only in the regime of strong magnetic quantization and 
weak degeneracy. Therefore it is convenient to write the 
generalized Saha equation using Eq. (Q) and describe 
the deviations from it by a separate factor A. For the 
ideal free energy of protons, we use Eq. (|35|), and the one 
for the atoms is given by Eq. (52). Thus the generalized 
Saha equation reads 



where 



[1 - eM-P^^cp)] Z^e^, (54) 



A^P^^,~H27^alX,n,)+p(^^-^] (55) 
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allows for deviations of the exact value of F-^^'' from that 
given by Eq. due to electron degeneracy and popula- 
tion of excited Landau levels. The distributions of atoms 
over the discrete quantum numbers and over the trans- 
verse pseudomomenta are given by Eqs. (|4^ ) and (48), 
respectively. 

The occupation probabilities can be presented as a 
product of two terms that arise from and -Fhs ■ 



(56) 



Hereafter, we exclude from our formulas by explicit 
use of the electroneutrality condition N, 
Coulomb factor reads 



Np. Then the 



In^ 



3/. 



c 



9 In r dlnrs 



(57) 



where fg{rs,T) = PFg/{2Np) is described in Sec. pIB 
In the Debye-Hiickel limit, w'^ is given by |p2|| 



Inujgn 



(58) 



The hard-sphere factor reads 



where Inw^o) = {d/dNp - d/dN^)[4:Ntotv] is the low- 
density limit of Inw^^ and rj is the packing fraction (18) 
Explicitly, 



4:7r 

ln«;(°)(ifi) = -— {(nn + np)lliK^) 



+ 3n^[lsAK^){P) + ll{K^){l)]} , (60) 
[Nl{{l')+3{l'){l))+N^Np{l') 



SNtotV 

+NuNu, ((?') + 3(P) /h. + 3(0 ll^ + /3 J 

+ NpNnA,+mljl,^], (61) 



where 



(62) 



The dissociation equilibrium is given by Eq. <\25\j, 
where Zu,2 is replaced by wn^ exp(xH9 ) and is mod- 
ified according to Eqs. ^) and (|5l|). From Eq. ( |2^ ) we 

obtain — (w^yw^^^heie 



jjs (1 -??/2)ln<^^- 15772 + V 



In w 



(0) 



47r 



+ {nn + np + 671h2)'h2] • 



(63) 



(64) 
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FIG. 5. Distribution of atomic occupation numbers at 
T = lO'* K for the magnetic field strengths B = 10^^ G 
and 10^"^ G (indicated). For each value of B, the distribu- 
tion over the quantum number s at i/ = and over !/ at s = 
is shown for two density values: p — 0.001 gcm""^ (hatched 
histograms) and 0.1 gcm~^ (shaded histograms). 



A solution of Eqs. (||)-(| 
ichiometric constraint 



!J), supplemented by the sto- 
i-p -r riB. + 2nH2 = '^0, yields the 
equilibrium abundances and the free energy value. The 
solution is sought by an iteration procedu re, in analogy 
with the zero-field case described in Sec. II C , and the 
EOS is obtained from Eq. (^). In the strongly quan- 
tizing magnetic field and in the nondegenerate regime, 
the EOS is a sum of three analytic terms: the ideal term 
Pid — ntotkBT + (4(t/3c) T'*, the contribution due to the 
Coulomb nonideality given by derivation of the fit de- 
[[II B , and the hard-sphere contribution 
r/y^ntotkBT. 



scribed in Sec 
Pm = 4r/ (1 - r;72)(l 



V. RESULTS 
A. Distribution of plasma particles 

1. Occupation numbers 

Figure ^ displays the distribution of the atoms over 
quantum states given by Eq. (^9|) at i? = 10^^ Q and 
10^^ G, T — 10^ K, and at two relatively low densities 
p = 0.001 gcm^'^ and 0.1 gcm^"^. The left panel shows 
the relative occupation numbers for the tightly bound 
states v = 0, for different quantum numbers s. The dis- 
tribution is broader for higher density. This apparently 
surprising feature is easily explained by the presence of 
the third quantum parameter K±, in addition to s and v. 
At low density most atoms reside in the states with large 
values of K± because of the large st atistica l weight of 
such states, which all have s = (Sec. IV A 1). At higher 
density, these strongly decentered states arc removed by 
the excluded-volume effects, and the distribution over s 
grows broader. Conversely, on the neighboring panel we 
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observe a narrower distribution over v at higher density, 
because the excluded- volume effects eliminate the hydro- 
genlike states. In the next section we shall see that ulti- 
mately, at still larger densities, only the ground centered 
state survives (s = i/ = 0, < Kc). 

The right two panels demonstrate the effect of increas- 
ing B to 10^^ G. Due to the larger binding energies, the 
distribution at p = 0.1 gcm~^ has become narrow, with 
more atoms concentrated in the ground state. However, 
at the lower density, the distribution over v has changed 
weakly, since the increase of binding energies is accom- 
panied by a decrease of the atomic size (and hence a 
decrease of the nonideality effects). 



2. Ionization equilibrium 

Figure || shows the ionization curves at three values 
of T for B ~ 10^^ G. The heavy solid lines represent the 
total fraction of atoms /h = n-H/no in all quantum states, 
calculated according to Eq. (p4h. Thin solid lines show 
the fraction foo of atoms in the ground state {s = v = 0, 
but any K±), and the dashed lines show the fraction of 
atoms in the centered states {K± < Kc, any s and v). For 
reference, triangles display the zero-field atomic fraction 
given by Eq. (pO|). 

We see that a strong magnetic field generally increases 
the neutral fraction. At low densities, the excited atoms 
contribute significantly. Since their effective size is pro- 
portional to K± , the integration (pO|) gives roughly Zg^ cx 



-2/3. 



therefore /oo decreases asymptotically as Uq 



1/3 



Because of the broadening of the ly distribution (roughly, 



max^" oc n, 



-l/6x 



g ), the low-density wing of the curve for 



the total neutral fraction has a slope /h cx: n^^, which 

1 /2 

is very moderate compared to /h oc tiq in the nonmag- 
netic case (triangles). 

The centered atoms, whose pseudomomentum is lim- 
ited from above by the critical value have a nearly 
density-independent IFF at low p. Therefore their frac- 
tion behaves as fcen oc no, and they disappear much 
faster at low p and especially at high T (compare the 
dashed lines in the upper and lower panels). 

At high densities, on the contrary, the decentered 
states become depleted due to the excluded-volume ef- 
fects, so that the dashed line in the figure merges with 
the solid one at p > 10 gcm^'^. At still higher densities 
p > 10^ gcm~'^, all excited states disappear, and only 
the state s = i' — survives. The pressure ionization 
proceeds around p ~ 10^ — 10'' gcm^'^. The excluded- 
volume and Coulomb nonidcal effects favor pressure ion- 
ization [both and in Eq. ( |5^ ) are less than 
unity], whereas finite electron degeneracy hampers it 
(A > 0). Because of the reduced atomic volume, the 
pressure ionization occurs at densities orders of magni- 
tude larger than for the zero-field case |^^. At T = 
10^'^ K, the molecular fraction becomes non-negligible 
at p ^ 10^ gcm^^. 
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FIG. 6. Ionization isotherms at B = 10^^ G and three 
values of T (indicated): total fraction of atoms /h ~ nu/no 
(heavy solid lines) and the fractions of ground-state atoms 
(thin solid lines), the centered atoms (short-dashed lines), and 
the optically identifiable (Inglis- Teller) atoms (long-dashed 
lines). The dot-dashed lines show the molecular fraction 
= 2nH2/wo, which is below the frame in the bottom pan- 
els. For comparison, /h in the zero-field case (triangles) and 
in the approximation of Lai and Salpeter (dotted lines) is also 
shown. 

Not all of the neutral atoms that contribute to the 
EOS may be identified spectroscopically. Because of their 
perturbation by plasma microfields, the atoms that do 
not satisfy the Inglis- Teller criterion form "optical con- 
tinuum." An approximate estimate of the fraction of 
atoms below the optical continuum is given by a gener- 
alization of the optical occupation probabilities Wk (Sec. 
II D) to the case of the strong magnetic field according 
to Eq. (14) of Ref. This "IT" fraction is shown 

by the long-dashed lines. Their rapid decrease indi- 
cates that the atomic spectral features disappear around 
p ~ 10 gcm^"', long before pressure ionization. 
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FIG. 7. The same as in Fig. | for B = lO" G. The vertical 
hne in the top panel separates the region of thermodynamic 
instability. 

The approximation of Lai and Salpeter ||20|| , also shown 
in the figure (dotted line), clearly underestimates the 
neutral fraction at low density and overestimates it at 
high density, especially at high temperature. At low den- 
sity, the discrepancy arises mainly from an underestima- 
tion of the decentered states because of an incorrect fit- 
ting formula to their binding energies. From comparison 
with the dashed line in Fig. ^, we see that the fraction 
of centered states can be estimated by the approxima- 
tion ||2§ at p < 0.1 gcm-3 and T < 10^ K. At higher 
T or larger p, the atomic abundance is overestimated in 
^ , pO| because of neglecting nonideal effects. Although 
the neutral fraction is very significant, it never dominates 
the plasma at the values of T and B shown in Fig. con- 
trary to the prediction of Ref. |2^]. (At T = lO^-^ K, the 
maximum is /h = 0.41 at p « 5 gcm~'^.) 

Figure |^ shows the ionization curves for a stronger 
field, B = 10^^ G. Under this condition, the neutral frac- 
tion still increases. At T = 10^'^ K (top panel), /h ex- 
ceeds ^ at p > 0.1 gcm~^, reaching the maximum of 85% 
at p « 10. Most atoms in this regime reside in the cen- 
tered ground state. On the other hand, at T = 10^'^ K 
and p ~ 10^—10"^ g cni^^^, the molecules are the dominant 
species; hence our present model may be not accurate in 



this p — T domain. 

A comparison with the result by Lai and Salpeter is 
not performed for T = 10^'^ K (the bottom panel of 
Fig . 0) because the approximations (3.11), (3.12) of Ref. 
pQI yield a negative IPF in this case. 

At T = 10^-^ K and p > 300 gcm'^, there ap- 
pears thermodynamic instability {dP/dp < 0) leading 
to a phase transition. The stability is recovered at 
p ^ 8000 gcm"'^, where the plasma is fully ionized. This 
phase transition is a complete analogue to the plasma 
phase transition (PPT) predicted in the zero-field case by 
several theoretical models [|4|j6^ but not yet confirmed 
in experiment. It is caused by a strong Coulomb attrac- 
tion between pressure-ionized plasma particles, which 
contributes negative pressure that cannot be compen- 
sated at low temperature until the degeneracy sets in. 
There is no confidence in the reality of the PPT be- 
cause of its model-dependence. In our case, an additional 
uncertainty is introduced by the simplified treatment of 
molecules. 
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FIG. 8. Dependence of the atomic fraction /h (solid lines) 
and the fraction of the optically identifiable atoms (dashed 
lines) on the magnetic field strength at two values of p (indi- 
cated in the figure), T = 10^-^ K (upper curves and 10^-^ K 
(lower ones). The atomic fraction at B = (triangles) and 
the approximation of Ref. (dotted lines) are shown for 
comparison. 

The _B-dependence of the atomic fraction at two values 
of T and two values of p is shown in Fig. ^. The total /h 
is drawn by solid lines, and the "optical" (Inglis- Teller) 
fraction by dashed lines. Triangles in the left panel show 
the total fraction of atoms at _B = (it is negligible at p = 
10 gcm~^ on the right panel). Dotted lines correspond 
to the approximation Q at T = 10^'^ K. 

It was found previously |l^,|l^ that the ionization de- 
gree decreases with growing B above ^ 10^^ G only at 
T < 5 X 10^ K but, in contrast to the present results. 
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increases at higher T. This behavior was attributed 
to two effects: decreasing phase space occupied by a 
plasma particle with growing B, which favors ionization, 
and increasing binding energy, which disfavors it. Our 
present result arises from the motional perturbations of 
the atoms, neglected in first, increasing B in- 

creases the effective mass nr^ and thus the statistical 
weight of the centered atoms, and second, at low densities 
the atomic IPF is further increased due to the decentered 
states. 



B. Equation of state 

Figure ^ presents four pressure isotherms obtained us- 
ing the free-energy model described in Sec. IV. For com- 
paris on, w e also show the fully ionized ideal-gas EOS 
(Sec. |III A| ) and the nonmagnetic EOS (Sec. The ver- 
tical line bounds the region p < ps- Let us first discuss 
the low-density regime p < 10 gcm~'^. At T > 10® K, 
all three EOS reduce to P = ngfesT. At lower temper- 
atures, the pressure deviates from this law because of 
the partial recombination of atoms. As discussed in the 
previous section, a strong magnetic field increases the 
neutral fraction; therefore the pressure is reduced more 
significantly compared to the B = case. 

In the intermediate-density range 10 gcm~'^ ^ p ^ ps, 
the differences among the three considered cases are most 
important. For B = 0, the plasma is fully ionized in this 
region, and the electrons become partially degenerate, 
making the EOS stiffer. In a strong magnetic field, the 
electron degeneracy is reduced (Sec. Ill); hence the ideal- 
gas EOS is softer, except for densities approaching pB, 
where the degeneracy sets in and pressure grows rapidly. 
Partial recombination and Coulomb nonideality lead to 
still further decrease of P. The pressure ionization dis- 
cussed above has two opposite effects on the pressure: the 
positive ideal-gas contribution of free electrons appear- 
ing in the course of the ionization and the positive non- 
ideal pressure of neutral species compete with the nega- 
tive Coulomb contribution. At low T, these effects may 
cause the thermodynamic instability mentioned above, 
which we observe on the isotherm T = 10^ K. The sec- 
ond lowest isotherm in the figure is slightly overcritical 
for this PPT. The dependence of the critical temperature 
Tc and density pc on B can be fitted by simple power laws 
Tc = 3 X 105B?2^9 K and pc = USBji^ gcm'^, where 
Bi2 = i?/(10^^ G). These fits provide an accuracy of a 
few percent in the considered range of the field strengths 
7 X 10" G < S < 3 X 1013 G. 

At higher density p ^ pb, excited Landau levels be- 
come populated due to the increase of the Fermi energy. 
Eventually, at p S> ps, the nonmagnetic EOS is recov- 
ered. 
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FIG. 



9. EOS of partially ionized atomic hydrogen at 
B = 10^^ G (solid lines) compared with the EOS of 
fully ionized ideal electron-proton plasma (dotted lines) and 
the EOS of partially ionized hydrogen at B = (dashed 
lines) . The temperature logarithms are (from top to bottom) 
\ogT[K\ — 6.5, 6.0, 5.5, and 5.0. The vertical line corresponds 
to pB , above which excited Landau levels become populated. 
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FIG. 10. Density exponent Xp ~ {din P/d\n p)t at 
T = 10*^ K without magnetic field (dashed line) and in strong 
magnetic fields of various indicated strengths (dot-dashed and 
solid lines). 

Figure ^ demonstrates the effects of the strong mag- 
netic field on the density exponent Xp — {d\nP/d\np)T- 
Although the pressure approaches the nonmagnetic value 
at p > pb, the effects of magnetic quantization remain 
quite prominent for the derivative Xpi shown by the 
curve B = 10^^ G in the figure. Consecutive population 
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of excited Landau levels causes the oscillations of Xp and 
other second derivatives of F around their nonmagnetic 
values. The regime where these oscillations are signifi- 
cant is called weakly quantizing 



1.5 



1 1 1 1 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 
logic ^ = 6 " 


-log.o'B = 




^-^^ ^^^^ 

X \^ 
\ ^^^^ - 

\ 

\ 

\ 




12.5 ' ' 


- - ' '^ v ^■ 






1 1 1 1 


12.0 

1 1 1 1 1 1 


\ - 
\ 

1 1 1 1 1 1 1 1 1 1 1 1 1 



-2-10 1 2 3 

logioP (g cm-3) 

FIG. 11. Normalized heat capacity at the same T and B 
values as in Fig. |l^. 

The effects of a strongly quantizing magnetic field on 
the reduced heat capacity Cy = k~^(dU jdT^y divided 



by the number of plasma particles, Aft( 



A^H + .^H2 1 a-re shown in Fig. In the nonmagnetic 
case (dashed line), the classical value Cy/A^tot = § is 
slightly exceeded at lower densities because of the ther- 
mal ionization of the atoms, and it is reduced to smaller 
values at higher densities because of electron degeneracy. 

In strong magnetic fields, the heat capacity is modi- 
fied due to several effects. In the low-density regime, Cy 
is reduced compared to the nonmagnetic value because 
of the quantization of the transverse motion of electrons 
and protons. The strongly quantized electrons have only 
one motional degree of freedom, so that their contribu- 
tion reduces to Cy^ — \^e- When protons are non- 
quantized and the plasma is fully ionized, this amounts 
to Cy /N tot = 1. In the general case, the contribution of 
free spinless protons would be 



(1) 

Vp 



2 smh{^ PTiuJcp) 



Nr. 



(65) 



which tends to ^ at jdhucp = 0.732 _Bi2/r6 ^ 1, where 
the protons are strongly quantized. The interaction of a 
magnetic field with proton spin, according to Eq. (36), 
yields 



^Vp - 



PgpfiU). 



cp 



A cosh{ jf3gpnuJcp) 



Nn 



(66) 



which vanishes in the limiting cases of huJcp ^ ksT and 
fiLOcp ^ ksT. In the latter case, Cy/Ntot would tend to i 
for the fully ionized plasma. In Fig. [l^, however, this does 
not happen because of the contribution of neutral atoms, 



which are subject to thermal ionization in this p-T-B 
domain. On the contrary, Cy increases with increasing 
B, since the neutral fraction becomes larger. The two 
humps visible on each magnetic isotherm correspond to 
the regions of the pressure destruction of the first excited 
atomic state s = 1, = and the ground state s — 
1^ = 0, respectively. In the latter case, Cy even exceeds 
the nonmagnetic value, because of the delayed onset of 
degeneracy. Only with density approaching ps, is the 
zero-field value of the heat capacity recovered. 

This illustrates the main effects of a strongly quantiz- 
ing magnetic field on a partially ionized hydrogen plasma. 
Other thermodynamic quantities, obtained within the 
framework of the present model, experience similar pro- 
found modifications. 



VI. CONCLUSIONS 

We have developed a thermodynamic model of the hy- 
drogen plasma in strong magnetic fields, making use of 
the available quantum-mechanical results for the fully 
ionized plasma and for the hydrogen bound species. Ap- 
plicability of the developed model is limited to the tem- 
peratures T, densities p, and magnetic field strengths B 
at which formation of molecules and other bound species 
more complex than the H atoms may be neglected. This 
condition holds, for instance, at B12 ^ 10 and T > 10^ K 
(any p) or at T > 10^ K and p < lO'' {Te/Bi2f gcm~^. 
Furthermore, although the theory presented in Sec. ^ 
is rather general, our numerical results in partially ion- 
ized regions are restricted to B12 > 0.7, because fitting 
formulas for quantum-mechanical characteristics of 
the atoms moving in magnetic fields have been derived 
under this condition. This restriction is fulfilled for a 
majority of neutron stars. For laboratory field strengths 
(at 7 <g; 1), perturbative methods may be sufficient. 

Calculations in the frames of our model show that 
the magnetic field effects strongly modify the thermo- 
dynamic functions and phase diagram of the plasma, in 
particular the partial ionization region. The abundance 
of atoms is significant in the considered domain of tem- 
peratures T ^ 10^ — 10^'^ K and magnetic field strengths 
B - 10^2 -IQis G at densities up to p - 102-10^ gcm^^^ 
contrary to the zero-field case. At relatively low densities 
(p < 1 — 100 gcm~^, depending on B and T), the decen- 
tered atomic states possessing a large constant dipole mo- 
ment are significantly populated. Since these values of p, 
T, and B are typical of atmospheres of isolated neutron 
stars, the physical effects discussed above are expected 
to affect the spectra. It has been shown [^5[^ that the 
presence of a nonionized component and, in particular, 
decentered atoms should produce observable absorption 
and thus necessitate a modification of previous fully ion- 
ized atmosphere models [|68[. Work in this direction is 
under way |69|| . 



15 



ACKNOWLEDGMENTS 

A.Y.P. is grateful to the theoretical astrophysics group 
at the Ecole Normale Superieure de Lyon for hospitality 
and financial support. A.Y.P. and Y.A.S. acknowledge 
useful discussions with J. Ventura and (at the initial stage 
of the work) with G. G. Pavlov and V. E. Zavlin. We 
thank M. Steinberg for sending us some of his results 
prior to publication. The work of A.Y.P. and Y.A.S. has 
been partially supported by INTAS Grant No. 96-542 and 
RFBR Grant No. 99-02-18099. 



V. Canute and J. Ventura, Fundam. Cosm. Phys. 2, 203 
(1977). 

D. G. Yakovlev and A. D. Kaminker, in The Equation 
of State in Astrophysics, edited by G. Chabrier and 

E. Schatzman (Cambridge University Press, Cambridge, 
England, 1994), p. 214. 

G. G. Pavlov, Yu. A. Shibanov, V. E. Zavlin, and R. D. 
Meyer, in The Lives of the Neutron Stars, edited by M. A. 
Alpar, U. Kiziloglu, and J. van Paradijs (Kluwer, Dor- 
drecht, 1995), p. 71. 

G. Chabrier, A. Y. Potekhin, and D. G. Yakovlev, Astro- 
phys. J. 477, L99 (1997); A. Y. Potekhin, G. Chabrier, 
and D. G. Yakovlev, Astron. Astrophys. 323, 415 (1997). 
D. Page, Yu. A. Shibanov, and V. E. Zavlin, in 
Rdntgenstrahlung from the Universe, MPE Report 
No. 263, edited by H. U. Zimmermann, J. Triimper, and 

H. Yorke (MPE, Garching, 1996), p. 173; D. Page, As- 
trophys. J. 479, L43 (1997). 

D. H. Constantinescu and P. Rehak, Nuovo Cimento B 
32, 177 (1976); D. H. Constantinescu and G. Moruzzi, 
Phys. Rev. D 18, 1820 (1978). 

I. Fushiki, E. H. Gudmundsson, and C. J. Pethick, As- 
trophys. J. 342, 958 (1989). 

A. M. Abrahams and S. L. Shapiro, Astrophys. J. 374, 
652 (1991). 

O. E. Rognvaldsson, I. Fushiki, E. H. Gudmundsson, 
C. J. Pethick, and J. Yngvason, Astrophys. J. 416, 276 
(1993); A. Thorolfsson, 6. E. Rognvaldsson, J. Yngva- 
son, and E. H. Gudmundsson, ibtd. 502, 847 (1998). 
Yu. N. Gnedin, G. G. Pavlov, and A. I. Tsygan, Zh. 
Eksp. Teor. Fiz. 66, 421 (1974) [Sov. Phys. JETP 39, 
301 (1974)]. 

V. K. Khersonskii, Astron. Zh. 64, 433 (1987) [Sov. As- 
tron. 31, 225 (1987)]. 

V. K. Khersonskii, Astron. Zh. 64, 1233 (1987) [Sov. As- 
tron. 31, 646 (1988)]. 

L. P. Gor'kov and L E. Dzyaloshinskii, Zh. Eksp. Teor. 
Fiz. 53, 717 (1967) [Sov. Phys. JETP 26, 449 (1968)]. 
L. A. Burkova, I. E. Dzyaloshinskii, G. P. Drukarev, and 

B. S. Monozon, Zh. Eksp. Teor. Fiz. 71, 526 (1976) [Sov. 
Phys. JETP 44, 276 (1976)]. 



[15] M. Vincke, M. Le Dourneuf, and D. Baye, J. Phys. B 25, 
2787 (1992). 

[16] G. G. Pavlov and P. Meszaros, Astrophys. J. 416, 752 
(1993). 

[17] A. Y. Potekhin, J. Phys. B 27, 1073 (1994). 

[18] J. Ventura, H. Herold, H. Ruder, and F. Geyer, Astron. 

Astrophys. 261, 235 (1992). 
[19] D. Lai and E. E. Salpeter, Phys. Rev. A 52, 2611 (1995). 
[20] D. Lai and E. E. Salpeter, Astrophys. J. 491, 270 (1997). 
[21] M. Steinberg, J. Ortner, and W. Ebeling, Phys. Rev. E 

58, 3806 (1998). 
[22] W. Ebeling, W. D. Kraeft, and D. Kremp, Theory of 

Bound States and Ionization Equilibrium of Plasmas and 

Solids (Akademie, Berlin, 1977). 
[23] W. Dappen, L. Anderson, and D. Mihalas, Astrophys. J. 

319, 195 (1987). 
[24] A. Y. Potekhin, J. Phys. B 31, 49 (1998). 
[25] W. Dappen, Rev. Mex. Astron. Astrofis. 23, 141 (1992). 
[26] G. Chabrier, in The Equation of State in Astrophysics, 

edited by G. Chabrier and E. Schatzman (Cambridge 

University Press, Cambridge, England, 1994), p. 287. 
[27] E. Fermi, Z. Phys. 26, 54 (1924). 

[28] D. G. Hummer and D. Mihalas, Astrophys. J. 331, 794 
(1988). 

[29] D. Mihalas, W. Dappen, and D. G. Hummer, Astro- 
phys. J. 331, 815 (1988). 

[30] W. Dappen, D. Mihalas, D. G. Hummer, and B. W. Mi- 
halas, Astrophys. J. 332, 261 (1988); D. Mihalas, D. G. 
Hummer, B. W. Mihalas, and W. Dappen, ibid. 350, 300 

(1990) . 

[31] D. Saumon, G. Chabrier, and H. M. Van Horn, Astro- 
phys. J. Suppl. Ser. 99, 713 (1995). 
[32] A. Y. Potekhin, Phys. Plasmas 3, 4156 (1996). 
[33] D. Saumon and G. Chabrier, Phys. Rev. A 44, 5122 

(1991) . 

[34] D. Saumon and G. Chabrier, Phys. Rev. A 46, 2084 

(1992) . 

[35] F. J. Rogers, Astrophys. J., 310, 723 (1986). 

[36] L. D. Landau and E. M. Lifshitz, Statistical Physics, 

Part I (Pergamon, Oxford, 1986). 
[37] H. M. Antia, Astrophys. J. Suppl. Ser. 84, 101 (1993). 
[38] M. Bans and J. P. Hansen, Phys. Rep. 59, 1 (1980). 
[39] S. Ichimaru, H. lyetomi, and S. Tanaka, Phys. Rep. 149, 

91 (1987). 

[40] G. Chabrier, J. Phys. (Paris) 51, 1607 (1990). 
[41] G. Chabrier and A. Y. Potekhin, Phys. Rev. E 58, 4941 
(1998). 

[42] H. DeWitt, W. Slattery, and G. Chabrier, Physica B 228, 
158 (1996). 

[43] J. Jung, M. S. Jhon, and F. H. Ree, J. Chem. Phys. 100, 

528 (1994); 102, 1349 (1995). 
[44] N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 

635 (1969). 

[45] Following [^, we estimate In ~ oqu \/7n'^ + 5/2 for an 
atom with the principal quantum number n. This yields 
doo = ao VS, which is 13% smaller than the SC ^] hard- 
core diameters. 

[46] In the partially degenerate regime, convergence of this 
procedure may be poor. In that case we impose a sta- 
bilization restriction that ui^ does not vary appreciably 



16 



from its value at the preceding iteration. 

[47] F. J. Rogers and C. A. Iglesias, Astrophys. J. Suppl. Ser. 
79, 507 (1992); C. A. Iglesias and F. J. Rogers, Astro- 
phys. J. 464, 943 (1996). We are grateful to F. J. Rogers 
for providing us with the monochromatic opacity tables. 

[48] V. E. Zavlin, G. G. Pavlov, and Yu. A. Shibanov, Astron. 
Astrophys. 315, 141 (1996). 

[49] F. J. Rogers, F. J. Swenson, and C. A. Iglesias, Astro- 
phys. J. 456, 902 (1996). 

[50] L. D. Landau and E. M. Lifshitz, Quantum Mechanics 
(Pergamon, Oxford, 1976). 

[51] A. Y. Potekhin and D. G. Yakovlev, Astron. Astrophys. 
314, 341 (1996). 

[52] F. Cornu, Europhys. Lett. 37. 591 (1997); Phys. Rev. E 
58, 5268 (1998); 58, 5293 (1998). 

[53] R. W. Danz and M. L. Glasser, Phys. Rev. B 4, 94 (1971). 

[54] F. Perrot and C. Dharma-wardana, Phys. Rev. A 30, 
2619 (1984). 

[55] A. Alastuey and B. Jancovici, Physica A 102, 327 (1980). 

[56] H. E. DeWitt, M. Schlangcs, A. Y. Sakakura, and 
W. D. Kraeft, Phys. Lett. A 197, 326 (1995); see also 
A. Alastuey and A. Perez, Phys. Rev. E 53, 5714 (1996). 

[57] E. E. Salpeter, Astrophys. J. 134, 669 (1961). 

[58] It was recently suggested, however, that this same 
large separation in the decentcrod states may help to 
create a positronium atom stable against the annihi- 
lation in the laboratory [J. Ackermann, J. Shertzer, 
and P. Schmelcher, Phys. Rev. Lett. 78, 199 (1997); 
J. Shertzer, ,J. Ackermann, and P. Schmelcher, Phys. Rev. 
A 58, 1129 (1998)]. 

[59] U. Willc, J. Phys. B 20, L417 (1986); Phys. Rev. A 38, 
3210 (1987); U. Kappes and P. Schmelcher, ibid. 53, 3869 

(1996) . 

[60] D. Neuhauser, S. E. Koonin, and K. Langankc, Phys. 
Rev. A 36, 4163 (1987); D. Lai, E. E. Salpeter, and S. L. 
Shapiro, ibid. 45, 4832 (1992). 

[61] D. Lai and E. E. Salpeter, Phys. Rev. A 53, 152 (1996). 

[62] Yu. P. Kravchenko and M. A. Liberman, Phys. Rev. A 
55, 2701 (1997). 

[63] T. Detmer, P. Schmelcher, F. K. Diakonos, and L. S. 
Cederbaum, Phys. Rev. A 56, 1825 (1997), and refer- 
ences therein; T. Detmer, P. Schmelcher, and L. S. Ceder- 
baum, tbid. 57, 1767 (1998). 

[64] M. A. Ruderman, Phys. Rev. Lett. 27, 1306 (1971). 

[65] G. G. Pavlov and A. Y. Potekhin, Astrophys. J. 450, 883 
(1995). 

[66] G. E. Norman and A. N. Starostin, High Temp. 6, 394 
(1968); W. Ebeling and W. Richcrt, Phys. Lett. 108A, 
80 (1985); G. A. Baker, Phys. Rev. E 56, 5216 (1997). 

[67] A. Y. Potekhin and G. G. Pavlov, Astrophys. J. 483, 414 

(1997) . 

[68] Yu. A. Shibanov, G. G. Pavlov, V. E. Zavlin, and J. Ven- 
tura, Astron. Astrophys. 266, 313 (1992). 

[69] Preliminary results have been presented by A. Y. 
Potekhin, Yu. A. Shibanov, and J. Ventura, in Neutron 
Stars and Pulsars, edited by N. Shibazaki, N. Kawai, 
S. Shibata, and T. Kifune (Universal Academy Press, 
Tokyo, 1998), p. 161. 



17 



